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ABSTRACT 

We present a study of outflow and feedback in the well-known Seyfert 2 galaxy 
Markarian 573 using high angular resolution long-slit spectrophotometry obtained with 
the Hubble Space Telescope Imaging Spectrograph (STIS). Through analysis of the kine- 
matics and ionization state of a biconical outflow region emanating from the nucleus, 
we find that the outflow does not significantly accelerate the surrounding host-galaxy 
interstellar gas and is too weak to be a strong ionization mechanism in the extended 
emission regions. Instead, the excitation of the extended regions is consistent with pho- 
toionization by the active nucleus. From energetics arguments we show that the nuclear 
outflow is slow and heavy and has a mechanical luminosity that is only ~1% of the 
estimated bolometric luminosity of the system. The energy in the outflow is able to 
mildly shape the gas in the extended regions but appears to be insufficient to unbind it, 
or even to plausibly disrupt star formation. These results are at odds with the picture 
of strong AGN feedback that has been invoked to explain certain aspects of galaxy 
evolution. 
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Introduction 



Recent studies of galaxy evolution have invoked a significant component of energetic feedback 
from an Active Galactic Nucleus (AGN) to attempt to explain many of the observed properties 
of galaxies. Specifically, AGN feedback has been implicated as the primary cause behind the 
differences between observed and theoretical galaxy luininosity functions, the observed color bi- 
modalitv of galaxies (jwhite fc Frenkll l99l': 'S Dringeil2nnnl : benson et all 1200.4 branato et ahlbooi 



Kauffmann et al 



2004 



Springel et al.l i2005al ) in which galaxies are divided into dis tinct popula- 



tions of blue star-forming galaxies and red dead galaxies (e.g. IStrateva et al.l l200ll ) with a pro- 
nounced deficit of galaxies with intermediate color, and the very tight correlations observed be- 
tween the masses of supermassive black holes an d their host gal axy spheroid's velocity dispersion 
(jPerrarese Merritt|[200ol : bebhardt et ahlboool ) and luminosity ^Marconi Huntlbood ). The rel- 
ative lack of galaxies with inter mediate star fo rmation rates suggests that star formation ceases 
abruptly rather than gradually (iBell et al.ll2004l ). while the tight correlations of black hole masses 
with the properties of their host bulges are taken as evidence that the two are strongly co-eval. In 
the various proposed feedback scenarios, the supermassive black holes grow primarily by accretion 
until they become sufficiently massive and energetic that thermal and/or radiative feedback from 
their ac tivity either heats the s urrounding inter stellar gas or mechanically pushes it out of the host 



galaxy (jCrenshaw et al 



2003 



raw m aterial from which to form new stars (jSpringel et al 



BegelmanI l2004l). This simultaneously deprives the host galaxy of 



20061 ) and starves the black hole, truncating its growth (jSilk &: Rees 199 



2005^ : iBower et al 



2001 



Fabian 



i 



Croton et al 



199 



J. 



Despite the considerable utility of AGN feedback, there is scant observational data to inform 
us about how or even if AGN feedback operates in actual galaxies. Arguably some of the most 
promising cases are instances of mechanical feedback, such as the outflows that have been de- 
tected as blue-shifted UV and X-ray intrinsi c absorbers in ~50% of Seyfert 1 galaxies and quasars 



( Crenshaw et al 



19991 . 



2003 



Dai et al.l 120081 ). While it is clear that nuclear outflows are common, 



it is still difficult to determine the mechanical luminosity (critical for determining if the outflow can 
unbind circumnuclear gas) and mass outflow rate to determine their impact. A good example of 
this difficulty is the work of iKrongold et al.l (j2007l ) on the time evolution of the ionization state of 
the X-ray absorbers relative to the X-ray continuum in the warm absorber of NGC 4051. From the 
absorption line variability they obtained good constraints on the density and location of the ab- 
sorbers from the black hole. Assuming a biconical geometry, they measured a low outflow velocity 
with respect to the escape speed from the black hole, and a corresponding low mass outflow rate 
with respect to the accretion rate of the black hole. From these values, they conclude that while the 
outflows might disr upt the hot ISM, they are not ca pable of ejecting large amounts of interstellar 
gas from the host. IWhittle et al.l (120021 . l2005l . l2009l ) reached a similar conclusion for the nearby 
Seyfert galaxy Markarian 78 using visible-wavelength and radio data to measure the kinematics and 
ionization state of extended emission-line regions surrounding the active nucleus. They concluded 
that the biconical outflow in Mrk 78 is weak, slow and heavy, and hence insufficiently energetic to 
significantly heat or unbind surrounding host galaxy gas. 
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The results for these two galaxies present an apparent problem for AGN feedback: in both cases 
the measured AGN outflows are too weak to thermally or mechanically disrupt star formation in 
the host galaxy. Nevertheless, these are just two objects, and analyses of other AGNs are critical for 
further understanding the role of AGN feedback in galaxy evolution. Seyfert galaxies are of special 
interest for detailed studies of AGN feedback as many are luminous enough to be accreting near their 
maximal Eddington rate, yet close enough to study at relatively high angular resolution (scales of 
10-100 pc for nearby examples) using the Hubble Space Telescope and ground-based adaptive optics. 

A very promising candidate for detailed study is the nearby Seyfert galaxy Markarian 573, 
which is well-known for its extended, richly structured circumnuclear emission-line reg ions. Mrk573 



features a prominent ioii ization bicone and bright arcs and knots of emission-line gas ([Ferruit et al 



199S 



Quillen et al.lll999l^ that are strongly aligned and interacting with a kilopars ec-scale low-power 



radio outflow (jPogge fc de Robertia Il993l : iFalcke et al.l Il998l : iFerruit et a. 



19991^. The assoc iation 



between the radio lobes and the emission regions was studied in detail by IFalcke et al.l (jl998l ). who 
proposed that the arcs result from gas c ooling after passing through a radio-induced radiative bow 
shock. In contrast, iQuillen et al.l (jl999l ) argued that the arcs are instead a morphological artifact 
of dust lanes in the galaxy being illuminated by the nuclear ionizing continuum. This is supported 
by visible and near-infrared HST images that show the emission-line arcs as extensions of larger 
circumnuclear dust lanes that are illuminated when they pass into the ionization co ne. Circumnu 



clear dust lanes like these are commonly found in active and inactive galaxies alike (jMartini et al 



20031') . and in the sp ecific case of Mrk 573 they appear to be shaped by the presence of a nuclear bar 
(jMartini et ahlboOll ). not unlike what is seen in galaxies with inner bars regardless of the presence 
of nuclear activity. 



Ferruit et al.l (|l999l ) examined several possibilities for the origin of the arcs and extended 



emission-line regions using a combination of HST images and ground-based, integral-field spec- 
troscopy. They modeled the excitation in the arcs as either shock features, linked to radio jets, or 
pre-existing structures photoionized by the nucleus. Using standard emission-line diagnostics, they 
showed that the inner arcs are excited by the central continuum source rather than by fast, pho- 
toionizing shocks, citing as evidence that they detected no signs of a strong kinematic interaction 
through a radio-induced shock. The excitation of the outer arcs was not completely explained by 
nuclear photoionization; these arcs require an external source of photons in addition to the nucleus 
to account for the excitation levels. However, their analysis suggested that both the inner and 
outer arcs show little evidence that they are shock excited. 

In this paper we present new long-slit spectrophotometry of Mrk 573 obtained with the Hubble 
Space Telescope Imaging Spectrograph (STIS). The high angular resolution of STIS, combined with 
good velocity resolution in the bright Ha+[N II] emission lines allows us to examine the excitation 
and kinematics in these regions in greater detail than possible in previous studies. Our analysis 
of Mrk 573 follows a threefold approach. First, we examine the kinematics of the circumnuclear 
emission-line regions to separate kinematically-disturbed gas from quiescent gas in the rotating disk 
of the galaxy (see sec. l3.ip . The existence of these two kinematical components was clear from the 
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earlier ground-based work, but the greater angular resolution provided by Hubble lets us pinpoint 
the regions of disturbed and undisturbed gas on a scale of 10s of parsecs. Next we measure the 
densities and temperatures in the emission-line regions using standard nebular diagnostic emission 
lines. We examine the excitation state of these regions, comparing them with the predictions 
of published photoionization and shock excitation models to clarify the nature of the ionization 
mechanism in the bright emission line areas (see sec. l3.2p . Finally, we u ndertake a qu a ntitat ive 
analysis of the outflow energetics of Mr k 573 using techniques developed by I Whittle et alJ (120091 ) to 
examine the degree of feedback on the host galaxy surroundings (see seed]). Throughout this paper 
we adopt a distance of 74Mpc for Mrk573 {vheiio = 5150 it 11km s~^, Hq = 70km s~^Mpc~^), 
giving us a projected linear scale of ~360 pc arcsec"^. For reference, this is approximately half the 
distance of Mrk78 (vheiio = 11137km s"^). 



2. Observations and Data Reduction 

2.1. Hubble Space Telescope Observations 

We acquired spectra of the nucleus of Mrk573 on UTC 2001 October 17 using the Hubble 
Space Telescope and STIS, with the 52x0.2 aperture. Our data were obtained in one observation 
period lasting two orbits with two spectral settings: the medium-dispersion G750M grating centered 
on the Ha emission line, ranging from AA6300 — 6850A and the low-dispersion G430L grating 
providing coverage from AA2900 — 5700A. These configurations provide FWHM spectral resolutions 
for extended sources of 2.2 and 10. DA respectively, and angular resolution of O'.'OSl pixel^^. The 
spacecraft roll angle was unconstrained in our observation planning as the data were acquired as 
part of a program to observe the nucleus proper, and an attempt to constrain the slit position angle 
to the apparent jet axis resulted in plan windows too short for practical scheduling. By pure good 
fortune, the slit position angle was —70.78° , which is within 16° of the radio-axis position angle 
of — 54°but still passes through the bright part of the inner arcs, well within the ionization cone 
and the main region of the nuclear outflow seen by others. Figure [T] shows the location of the slit 
relative to the emi ssion-line regions on a c ontrast-enhanced "structure map" created from an HST 



F606W image (see lPogge Martinill2002l ). 



The G750M spectra were acquired as a sequence of three exposures with target exposure times 
of 1080, 1080, and 840s. The G430L observations consisted of two exposures of 805 and 840s 
duration. To avoid problems with hot pixels on the STIS CCD, the observations were dithered by 
±0'.'25 (5 pixels) along the slit relative to the first spectrum. Wavelength calibration lamp spectra 
were taken during Earth occultation. 
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2.2. Data Reduction 

Our data were acquired after tire primary (Side-1) STIS electronics failed on UTC 2001 May 16. 
Since the Side-2 electronics did not provide closed-loop temperature control of the STIS CCD, the 
CCD temperature varies with t he ambient temperature o f the spacecraft when the thermoelectric 



cooler is run at a constant rate (jKim Quijano et alJ 120031 ) . The result is that the dark rate varies 



with temperature, and the standard dark calibration images used by the calstis data reduction 
pipeline are often a poor match to the actual dark rates in the hot pixels. To correct for this, we 
broke out of the calstis pipeline after the BASIC2D step and corrected the dark pixels individually 
before re-inserting them into the pipeline for the final wavelength calibration, flux calibration, and 
geometric rectification steps. In brief, a custom program scanned the Data Quality Frame (DQF) 
for each individual image and created a list of hot pixels. This list was used to fix all pixels in 
the science frame that had either negative data values or that deviated from the median within 
a 7x 1-pixel box centered on that pixel by more than 2a". For this last step, the long-axis of the 
analysis box is oriented along the dispersion axis of the detector. Once this hot-pixel surgery was 
completed for the three G750M observations, the second and third images in the set were then 
aligned with the first and co-added using the calstis ocrreject task to detect and remove cosmic 
rays. Remaining cosmic rays missed by this step were removed by hand using the median filter 
routine TVZAP in XVista. The final combined, cleaned, spectral image was then fed back into the 
calstis pipeline to produce the final calibrated long-slit spectrum. An analogous process was used 
for each of the two G430L spectra to create the final G430L long-slit spectrum, but because only 
two images were available, more hand cleaning of residual cosmic rays was required. This process 
yielded superi or results by comp arison with the default reduction pipeline. An analogous procedure 



is described in iRice et al.l (j2006l ). 



An additional consequence of using the STIS Side-2 electronics is that our long-slit spectra are 
affected by low-amplitude (2-3 ADU peak-to-peak) fixed-pattern noise with a ~3-pixel horizontal 
scale that takes the form of a regular herring-bone noise pattern when viewed at high contrast. 
Attempts to remove this component using fast fourier-transform filtering only made things worse. 
Its primary effect is to increase the effective readout noise of the device by about 1 e~ pix~^ for 
the gain=l mode of our data. When extracting faint outer regions of the long-slit spectra, we 
took particular care to be wary of this fixed-pattern component in interpreting our spectra. The 
final, calibra.ted lon g-slit spectra are shown in Figure[2l We have adopted the nomenclature of 



Ferruit et al.l (119991 ) for the specific emission-line regions we examine below. 



2.3. Spectral Extraction and Line Measurements 

The extended emission-line region of Mrk573 consists of a series of arcs and knots located 
northwest and southeast of the nucleus along the general direction of the radio outflow. Our STIS 
slit intersects the arcs but just misses a group of bright knots located along the main radio axis 
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as shown in Figured]). We shal l focus our attentio n on four specific regions: the SEl arc, SE2 arc 



the nucleus, and the NWl arc ([Ferruit et al.lll999l ). The arcs at regions SEl and NWl are located 



at a projected angular distance of ~2" on either side of the nucleus, and the SE2 arc is ~3" from 
the nucleus. There is another arc in the Northwest region that is roughly symmetric with the SE2 
arc, but the emission lines from this region are too faint and diffuse for detailed analysis with our 
spectra. 

We extracted two ID spectral data sets from the 2D spectra. The first set is integrated spectra 
of specific regions of interest for measuring the density, temperature, and excitation state of the 
line-emitting gas. For these we first created slit intensity profiles for the Ha+[N II] emission lines 
in the G750M spectrum and for [O III] A5007A in the G430L spectrum, and for adjacent line- 
free continuum regions bracketing these lines. The derived mean continuum profiles were then 
subtracted from the line profiles to produce pure emission-line profiles for each line. These profiles 
were then used in conjunction with the 2D spectra of the nuclear regions (see Figure[2]) to precisely 
locate the windows for extracting ID spectra of the regions we wished to examine. 

The second set of spectra were extracted from the G750M data for making kinematic mea- 
surements - line centroids and velocity widths - using the Ha and [N II] emission lines. For these 
we extracted ID spectra in contiguous increments along the slit moving radially outwards from the 
nucleus. In bright regions we extracted spectra from single-pixel apertures (0'.'2x0'.'05), and as the 
regions became fainter at larger radii, we increased the width of the extraction window up to a limit 
of 4 pixels wide (0'.'2x0'.'2) to improve the signal-to-noise ratio in the lines. This gives us nearly 
continuous spatial sampling along the slit for measuring line-of-sight velocities and emission-line 
velocity dispersions. 

After identifying individual emission lines, we measured the line parameters with the LINER 
interactive line analysis program developed at Ohio State. For each spectral line, we define a 
local continuum by averaging over adjacent line- free regions and then fit the lines with a single- 
or multi-component Gaussian to resolve blends or multiple velocity components if present. We 
derive the line centroid. Full- Width at Half Maximum (FWHM), and integrated line intensity for 
each emission- line component present. When fitting doublet lines of a single ionic species (e.g., 
[N II] AA6548,83A and [S II] AA6716,3lA) we fit both lines together, imposing physically-motivated 
constraints on the profile parameters to improve the fidelity of the fits. For the [S II] doublet, 
we constrain the relative line centroids in velocity space and match them in line width, but leave 
the relative intensities (which are correlated with density) unconstrained. For the [N II] doublet 
lines, which arise out of the same upper excited level, we constrain their relative line centroids and 
widths and further constrain their relative line strengths as dictated by the ratio of their radiative 
transition probabilities. This greatly improves the quality of the fits when deblending these lines 
from the Ha emission line that lies between them. Some of the emission lines appear double- 
peaked, particularly in the kinematic spectral extractions. In these cases we model the two peaks 
as Gaussians and apply similar rules for setting constraints on the physically related line pairs. 
Finally, if the line of interest is very faint but clearly detectable and unblended, we extract the line 
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intensity by direct integration without attempting to fit a line profile. The RMS uncertainty in the 
continuum fits and a simple model for the pixel-to-pixel signal to noise ratio is used to estimate 
the overall quality of the fit and the uncertainties on the best-fit line parameters. 

We make additional checks of our fit-derived line parameters in a number of ways. The total 
line intensities are compared to a numerical integration of the blend itself without fitting to a 
particular line shape. Line FWHMs from the fits are compared to direct measurements of the 
FWHM using profile tracing algorithms that assume no underlying profile shape (this works best 
for the brightest emission lines with the least blending and highest signal-to-noise ratio). Finally, 
the Gaussian fit derived line centroids for well-isolated lines are compared to intensity-weighted 
centroids derived from the central 5 pixels about the peak. All of these give us further estimates of 
the relative uncertainties, and help alert us to possible systematics that would be missed by simply 
examining the RMS residuals. 

The measured fluxes for various emission lines are listed in Table[TJ In general, line centroids 
and FWHM are measured to itO.lA or ~ 20 km s"^ for the Hq and [N II] emission lines, with some 
degradation in fainter regions. Measurement uncertainties for total line fluxes are typically ±10% 
in the brighter lines, increasing to ±20% in the faintest regions measured. 

3. Spectral Analysis 

We pursue two different analyses of our spectra. The first uses the kinematic spectra (sec. 
13. 1|) to derive outflow speeds and separate outflowing gas from ambient ISM gas in orderly rotation 
in the disk. Combining the outflow kinematics with the emission-line geometry seen in the direct 
images provides us with constraints on the conflguration of the outflow, allowing us to deproject 
the observed radial velocities and thus determine the actual physical outflow velocities. 

The second line of analysis uses emission line diagnostics from the spectrophotometric set of 
spectra (sec. 13. 2p to estimate the density, temperature, and excitation state of the gas in these 
regions. These are compared to published models of shock- and photo-ionized gas to determine the 
excitation mechanism in these regions. 

3.1. Emission Line Kinematics 

Figures [3] and H] show the radial velocities along our long slit for the Ha A6563A and [O III] 
A5007A emission line^. A number of distinct kinematic systems are present. The largest-scale 
pattern is a rotating disk extending out as far as we can trace ionized gas (roughly ±4"), with the 



^Radial velocity measurements for H/? A4861, [N II] AA6548, 6583, and [O III] AA4363, 4959A lines show consistent 
results and are available upon request. 
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NW side receding and the SE side approaching us and an amphtude of roughly ±125 km s~^. In 
the inner 1" on either side of the nucleus the kinematics are dominated by a coherent system of 
double- valued radial velocities (separated by nearly 300 km s~^) with a red/blue splitting pattern 
characteristic of a biconical outflow (e.g. ICrenshaw et al.l l2000l : iDas et al.ll2005l ). Finally, in the 
inner 0'.'5 around the nucleus we see various additional Ha velocity components with no obvious 
pattern associated with the inner parts of the outflow region. The high-velocity components of the 
biconical outflow pattern and the large-scale rotation pattern are also visible in the [O III] A5007A 
line-of-sight velocities (FigureH]) despite the lower dispersion of the G430L grating. 



In the biconical outflow models of ICrenshaw et al.l (|200d ) and I Das et al.l (|2005l ) emission from 
the front of the cone (the surface closest to the observer) is blueshifted and the back of the cone 
is redshifted. When the bicone axis is at high inclination angles (the cone axis near the line of 
sight) there is a strong blue-to-red asymmetry in the radial velocities, whil e at low inclinat ions (the 
axis near the sky plane), the radial velocity pattern is more symmetric (IDas et al.ll2005l ). In the 
context of these models, our radial velocity maps indicate that the biconical outflow in Mrk573 has 
a relatively low inclination with the NW cone directed toward us. Using the v elocities from our Ha 
emiss ion line, the observed ionization cone opening angle of 9c = 45° it 10 (jWilson &: Tsvetanov 
19941 ). and assuming a constan t velocity, yf,oi m aterial in a shell along the cone and a slit roughly 
aligned with the cone axis, the lDas et al.l (|2005l ) outflow model gives projected velocities of 



VRed = VfSm{Oc + ' 
VBlue = Vf Sm{ec - • 



200km/s 
lOOkm/s 



(1) 
(2) 



for the red and blue sides of the cone, respectively. The cone opening angle {6c) is related to these 
velocites and the cone inclination angle ((/>) by 



tan{9c) 



f VRed + VBlue 
^ VRed — VBlue 



tan(0) ~ 3tan(0) 



(3) 



which gives a bicone inclination angle of = 18° it 4° and Vf 
in Mrk573. We model this outflow geometry in Figure[5j 



220 km s ^ for the nuclear outflow 



Between about 1-1'.' 5 from the nucleus there is very little emission line gas, and so we have no 
measurements in these regions. This region between the nucleus and the two close emitting arcs 
likely represents the nearly evacuated "bubbles" resulting from the biconical outflow. Beyond about 
1'.'5 on either side of the nucleus a coherent pattern of velocities consistent with a rotating disk of 
relatively undisturbed gas extends out to the limits of our long-slit data. This sense of rotation 
(SE approaching, NW receding) is opposite the direction of the biconical outflow which helps to 
distinguish these components. At the location of the bright arcs (SEl and NWl), the kinematics 
are disturbed: the radial velocities split by ~100km s^^relative to the rotating disk pattern, and 
the velocity dispersion jumps dramatically in the SEl arc to nearly 300 km s~^. This can be seen 
in the ID radial velocity profile (Figure[3l middle panel) and in the 2D spectrum (Figure[2]) where 
the Ha, [N II] and [S II] emission lines are all suddenly broader. The effect is less obvious in the 
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fainter NWl arc, although in the 2D spectrum th ere is faint bu t notic eable broadening of the Ha 
and [N II]A6583A lines. These are the regions that iFerruit et alj (|l999l ) identified with the working 
surface between the radio outflow and ambient gas in the galaxy. The SE2 arc appears to have 
similar kinematic disturbance as SEl and NWl in the [O III] velocity map; however, it is not 
consistently double- valued in Ha, like the other two arcs (see Figure[3lll|). This implies that SE2 
is not significantly kinematically disturbed as SEl and NWl are. Throughout this analysis, we 
include calculations for SE2, but our primary focus is on the effect of an outflow as manifested in 
SEl and NWl. 



3.2. Spectrophotometry 

The spectra of the nucleus and bright, off- nuclear regions are remarkably rich in emission lines, 
as can be seen in Figure[6] where we plot the spectra of the nucleus and the brightest off- nuclear re- 
gions. Ofparticular interest are the lines of Ha, H/?, [N II]AA6548,6583A, [O III]AA4363,4959,5007A, 
and [S Il]AA6716,3lA which are used to estimate internal extinction, densities, and temperatures 
in the gas, and provide diagnostics of shock versus photoionization. In addition, we can observe a 
number of fainter high-excitation lines of interest, including [Fe X]A6375A and the [Ne V]A3426A 
line. 



3.2.1. Densities and Temperatures 



The densities of the emission -line gas are esti i nated using the [S II] AA6716,3lA doublet ratio 
and the theoretical calculations of lCai &: PradhanI (|l993l ). The gas temperatures are ca lculated us- 
i ng th e [O III] AA4363,4959, and 5007A emission lines following the method described bv lOsterbrock 
(|l989l ). Because the [O III] emission-hne ratio 1(4959) + /(5007)//(4363) is reddening-sensitive, 
we estimate the effects of extinction by assuming a simple screen geometry and using the H I 
Balmer decrement assuming Case-B recombination. The small differences between the observed 
and emitted [O III] line ratios in Table[2] indicate that reddening in these regions is negligible. 

Because the [S II] doublet density method is weakl y temperature dep endent and the [O III] 
temperature estimate is weakly density dependent (see lOsterbrockl Il989l ) , we have adopted an 



iterative approach, whereby we first estimate the [O III] temperature in the low-density limit, 
evaluate the [S II] density at that temperature, and then use that density to re-evaluate the [O III] 
temperature, iterating until the solutions converge to within the measurement uncertainties. The 
densities we find range from 760-2500 cm~^, and the temperatures range between 10000-14500 K. 
The resulting temperature and density estimates for each region and their uncertainties are given 
in TableEl 
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3.2.2. Emission Measures 

The arcs are high-surface brightness emission- hne features, but the "bubbles" between SE 1, 
NW 1, and the nucleus (B^^; and BjwK respectively) are very low surface brightness. As the bubbles 
have only very weak emission lines, we cannot estimate the temperature and density using standard 
metal-line diagnostics. Instead, we make an order-of-magnitude estimate of the mean density (rig) 
using the emission measure [EM) derived from the surface brightness of H/3 {I (HP)), an estimate 
of the line-of-sight thickness of the bubbles, D, and the volume filling factor /: 

£;M= /^^^f^f^ = / n,n,ds^{nl)Df (4) 
hi^Hpaj^'p (T) Jios 

We assume that the temperature is roughly constant, and that the bubbles are roughly spherical 
(i.e., as thick as they are wide) and have unit filling factor. A first estimate assumes a temperature 
of 15000 K, slightly higher than the temperature of the nuclear region, to evaluate the effective 
recombination coefficient, giving a mean density of ~lcm~^ in the bubbles (see Table[3]). At this 
temperature and density the gas in the bubbles is out of thermal pressure balance with the gas in 
the emission-line arcs just outside them by a factor of ~500-1500. The bubbles may instead be more 
like supernova bubbles with internal temperatures of ~ 10^ K; adopting T = 10^ K gives emission- 
measure derived mean densities an order of magnitude larger (see Table[3]), bringing the bubbles 
into rough thermal pressure balance. The latter temperature is consistent with thermalization of 
bulk motions with velocities as observed, i.e. ~ 200 — 300 km s~^. 



3.2.3. Nuclear Emission 

While the bubble regions have only very weak emission lines, the nuclear emission spectrum has 
many high ionization lines that are undetected in other emission-line regions, notably [Fe X]A6375A, 
[Ar V] A6435A, and [Ar IV] AA4711, 4740A ( see Figure[6|) . The presence of these lines indicates that 
the photons emitted from the nucleus have a sufficiently hard spectrum to substantially ionize the 
surrounding material. Some of these lines appear to be double peaked; in particular [Fe X] A6375A 
has a broad, flat top. This feature is best fit by two blended Gaussian lines, at A6483 and A6486A, 
respectively, corresponding to line-of-sight velocities of —46 and +92 km s~^. This line splitting 
likely represents highly excited gas at the base of the biconical outflow. 



3.2.4. Emission Line Diagnostics 



Shock-ionized gas will have different relative emission line strengths than gas photoionized by 
the active nucleus. Collisionally excited UV resonance lines provide the best discriminants between 
shock and photoioniza tion due to the higher temperatures expected in the shocked gas (jPopita 



2OO2I : lAllen et al.lll998l ). Unfortunately, we do not have near-UV spectra of Mrk573 and therefore 
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we must use visible- wavelength indicators. 



To evaluate the importance of shock excitation for the arcs, we compare the fluxes of our 
observed optical emission lines (Table[T|) to the predictions from models for ionizing shocks with 
precursors, photons trav eling ahead o f the shock front, and nuclear photoionization calculated by 
Whittle et alJ (120051 ) and I Allen et alJ ([200a). The shock models cover a range of shock speeds {vs = 
200- 1000 km s-i) and magn etic field strengths, B / \/n = 0— 10//G cm^/^ . We adopt two illustrative 
photoionization models from [Whittle et al.l (|2005l ). The first is a power-law photoionization model 
for optically thick gas with a density of 10^ cm~^ and a nuclear continuum with a spectral index 
of Q = 1.4 (where fu oc i/~"). This model covers photon energies from 50keV to 100 /xm with an 
ionization parameter range of logU = —4.0 to 0.5. The second is the ^Vdust" model that is similar 
to the first but includes a realistic dust component that accommodat es changes in the in ternal 
structure of dusty clouds in response to the incident radiation field from lCroves et al~ ( 2004 ). This 
second model is for a density of lO^cm"^, continuum spectral index of o; = 1.4, and a metallicity 
twice the solar abundance. 

Figure [7] shows four emission- line diagnostic diagrams that provide the cleanest discrimination 
between shock and central photoionization for the lines available in our spectra. The upper two 
panels show two diagnostics using the [O I] A6300A hue: [O III]/H/? vs. [O IIl]/[0 I] and [O III]/H/? 
vs. [O l]/Ha. In both cases the data are in best agreement with the simple, Uia photoionization 
model and appear to disfavor the shock sequences and to a lesser extent the Udust model. However, 
these results should be viewed with some caution as the [O I] line is difficult to model: in photoion- 
ization models the [O I] is formed in a partially ionized zone heated by X-ray photons, while shock 
models are known to overpredict the [O I] line strength by as much as an order of magnitude (e.g. 
Groves et al1l2004l ^. 



The lower left panel of Figure [7] shows the [O IIl]/H/3 vs. [S II] /Ha diagnostic in which the 
photoionization and shock models are essentially orthogonal to one another. Here the nucleus and 
SE2 region are in better agreement with the photoionization model, while the NWl and SEl arc 
regions are ambiguous between photoionization and a relatively fast shock {v > 500 km s~^) with 
a strong magnetic field {B/^/n = 10//Gcm^/^), while they are in good agreement with the Udust 
model. In the last panel we show the [Ne V]/[Ne III] vs. [Ne IIl]/[0 II] diagnostic. These line ratio 
measurements are similar to the previous panel, in particular the agreement between the NWl and 
SEl measurements and the shock models, while the nucleus and SE2 region do not agree well with 
any of t he models. One ca veat with this diagram is that [Ne III] may be overpredicted by the 
models dwhittle et aklbood ). 



While none of these line ratios are ideal discriminants between shock and photoionization, 
we nevertheless see consistent trends in our analysis that favor photoionization. One interesting 
aspect of these data is that the SEl and NWl arc regions, which trace the most morphologically- 
obvious bowshock structures, are more consistent with shock models than either the nucleus or the 
SE2 region. This also supports our assertion that SE2 is not strongly interacting with a biconical 
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outflow. Additional far UV data could prove useful to further investigate the strengths of any 
shocks in these regions; however, based on these and previous arguments we conclude that the 
emission line ratios in all four regions are more consistent with photoionization than with shocks. 
While they themselves are not a clear cut indicator of photoionization, they do support the findings 
from kinematics and analysis of energetics. 



3.3. Discussion: Ionization Mechanisms 

Our analysis thus far adds additional detail to the picture for the inner regions of Mrk573 
developed by previous work: a hot biconical outfiow emerging from the nucleus that is running 
into the host galaxy at the location of the bright emission-line arcs. Our kinematic data show that 
this outflow is mildly inclined with respect to the host galaxy disk (18 ±4°), and there is a modest 
kinematic disturbance at the location of the arcs. The inner part of the outflow between the dense 
base of the bicone and the arcs has low density (~ lcm~^) 



The nature of the arcs from previous work has been equivocal. On the one hand, iFerruit et al 



(|l999l ) inconclusively modeled these same emission regions as either pre-existing structures or shock 



features, while on the other H ST images show dust lanes that extend beyond the ionization cones 



which led iQuillen et al.l (|l999l ) to suggest that the arcs of Mrk 573 do not result from bow shocks. 
Our analysis of the densities of the different emission line regions reveal that the arcs and nucleus 
have similar densities, implying that the arcs were not formed by a fast shock. A fast shock would 
have slammed into the surrounding material, making it more dense , which we do not observe. The 



density values support the morphological argument of lOuillen et al.l (jl999l ) tha t the central photon s 



are ionizing pre-existing dust lanes. In turn, this supports the suggestion of iFerruit et al.l (119991 ) 



that any shock from the central AGN is too weak and slow to photoionize the surrounding material. 

Analysis of the evacuated regions, the bubbles, compared to the emission line regions also 
supports central photoionization of dust lanes. If these bubbles were rapidly expanding and com- 
pressing the surrounding material, they would shock ionize the arcs. In this scenario we expect 
a higher pressure inside the bubbles than the arcs to cause the expansion. Our calculation of 
the emission measures (see sec. 13. 2. 2]) reveals that the bubbles have ~500-1500 times less thermal 
pressure than the arcs, assuming a temperature of 15000 K. When we assume a temperature of 
10^ K, motivated by plausible shock velocities, the thermal pressure in these regions is still only 
~0.1-1 times the amount in the emission line regions (see TableE]). For the bubbles to be rapidly 
expanding into the surrounding material, bulldozing it, they need to have higher pressure than 
their surroundings. Over a large range of temperatures, the pressure in the bubbles is simply not 
large enough for the amount of expansion that would result in arc creation and shock ionization. 
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4. Properties of the Nuclear Outflow 



Our spectrophotometric analysis of Mrk 573 indicates that the emission line regions are likely 
photoionized by the nuclear source. If there is a shock present, it is not strong enoug h to form 
the material into arcs and ionize it, agreeing with the conclusion of lFerruit et al.l (jl999l ). We now 
complement these analyses with a detailed investigation of the emission-line region and outflow 
energetics, quantifying the outflow's interaction with the host galaxy ISM. Specifically, we estimate 
the energies and pressures of the emission-line regions (section 14. ip and relate them to values 
calc ulated for the outflow (section 14. 2p . These analyses broadly follow the techniques applied by 
the IWhittle et al.l (|2009l ) study of Mrk 78, and we shall adopt similar notation. By comparing 
the values of the different energy components of the outflow and emission line regions, we can 
confirm that the emission line regions are photoionized. The various energies also indicate the 
importance of mechanical feedback, and the overall strength of the outflow present (section 14.31) 



We al so compare our results to those of the analogous emission-line regions studied bv I Whittle et al 



(|2009l ). and demonstrate that Mrk573's emission-line regions are very similar to what is observed 
in Mrk 78. 



4.1. Emission Regions: Masses, Energies and Pressures 

To understand the relationship between ionized and neutral gas in the emission regions, we esti- 
mate the values of the various energy components. This is critical for determining what contributes 
to the kinematics, ionization, and structure, in addition to the relative importance of mechanical 
feedback on the host. We are able to calculate the energy of the various components by building 
off of determinations of the flux, mass and lifetimes of the emission line regions. 

We assume that the portion of a given emission-line region within the slit is representative of 
the larger region as a whole, and so multiply the flux within the slit by the ratio of the total region 
area {6ex x Gey) to the slit region area to approximate the total flux. The slit regions are the width 
of the slit, 0.2" multiplied by the size in arcseconds each region covers in the pure emission-line 
profiles (see sec. 12. 31 Table[TJ. The physical sizes of the regio ns that we us e are l isted in TableSl 



Our region sizes and fluxes are smaller than those measured bv lFerruit et al.l (119991) becaus e we can 



better distinguish the sharp boundaries with our improved resolution; iFerruit et al.l (|l999l ) had an 



angular resolution of ~0'.'35 FWHM compared to our STIS resolution of G/.'l FWHM. 

The ionized gas mass in each emission- line region (Mem) may be estimated from the H/3 
emission-line luminosity and the electron density Ue derived from the [S II] doublet ratio from 
recombination theory: 

Mem ^ , ^^^""^ « 8 X 10^^ (^) VH^(nem,3)-' Mq (5) 



In this and subsequent equations, the numerical subscript appended to a variable indicates a log- 
arithmic scaling, e.g. riems represents density in units of 10^ cm~^. The interaction age of the 
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region, tem, is estimated as the crossing time, tcross, the amount of time required for a compact 
parcel of gas expanding outwards at Vem to traverse the distance from the nucleus to the arcs, i.e. 

tern ~ Across « 47100 ( ^]6ex{Ve.m,2)~^ tan (/) yr (6) 



■Ho 

where 9ex is the angular separation between the nucleus and arcs in units of arcseconds. The 
emission- line masses for NWl and SE2 are both ~ 10000 it 2200 Mq, while SEl is less massive, 
~ 7000 ± 1500 Mq. SEl has an interaction age of 3.53 ± 1.3 Myr, and NWl is 2.24 ± 0.7 Myr old. 
SE2 has the largest interaction age, tem ~ 6.90 ib 2.5 Myr. This is the expected age progression, in 
that SE2 lies further from the nucleus than SEl and NWl. However, we must note that SE2 does 
not appear to be strongly interacting with the biconical outflow, based on its lack of kinematic 
disturbance in Ha (see FigureE]). We use the interaction age of SE2 to determine other properties 
of the emission-line region to better understand the basic energetics of the region, but it is merely 
an approxim ation. The age deter minations for all of the regions are consistent with what is seen 



in Mrk78 bv lWhittle et al.l (|2009l ) from emission-line structures of comparable scale (0.4-8.2 Myr). 



The total energy in the emission-line regions, E^m^ may be estimated from the product of the 
interaction age and the luminosity in the [O III]A5007A emission line which measures the nebular 
cooling power: 



Eem = 1.2 X 10^1 (-^j>5007 X W (7) 

where -F5007 is the [O III] line flux in cgs units. Each of the regions has total emission energy on 
the order of 10^^ ergs. 

There are three energy components of interest in the emission line regions associated with their 
mechanical and thermal energy content. The translational kinetic energy, Ekin.t, of the outflowing 
gas may be estimated from the mass of the region and its transverse velocity, Vem, 

Ekin,t ~ ^MemV^^COSeC^^ « 10^^ Memfi {^em, 2)^ COSec'^ (j) (8) 

where we use the estimated bicone inclination angle (p to deproject the observed radial velocity. 
The internal (turbulent) kinetic energy, E^i^^i, depends on the emission-line mass and velocity 
full- width, Wem- 

Ekin,i ~ ^MemiWem/2.35)^ ^ 1.8 X W'^'' MemfiiWem,2? (9) 



2 



Whittle et al.l (120091) used the width of the [O III]A5007A emission line, but because we have only 
low-dispersion spectra in which these lines are either unresolved or marginally resolved at their 
broadest, we instead use the widths of the Ha and [N II] emission lines. Finally the thermal energy, 
Eth, is 

Eth ~ ^NkT ^A.lX l{f>^MemfiTe,i (10) 

where we use the electron temperatures derived from the [O III] emission lines. All of these energies 
are listed in TableUl In general, the translational kinetic energy in the regions is about 100 times 
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larger than the thermal energy. This is consistent with our earlier suggestion that the gas is not 
significantly shocked, or if shocked it has had enough time to cool to the ~ 10^ K temperatures 
observed. This implies that the acceleration mechanism must be fairly slow and gentle (the typical 
gas speeds are ~ 100 - 200km s"^ in these regions, about an order of magnitude larger than the 
sound speed for 10^ K gas). The internal kinetic energy is roughly 10 times larger than the thermal 
energy, suggesting that this motion is not supersonic turbulence but instead is bulk motion of 
the gas along the line of sight that is unresolved at our scales, for example an expansion of order 
iblOOkm s~^ along the line of sight at the point of contact between the nuclear outflow and the 
emission-line gas in the SEl and NWl arcs. 

Our spectrophotometric analysis of Mrk573 suggests that the emission-line regions are pho- 
toionized by the active nucleus rather than ionized by fast shocks. We investigate two sources 
of energy input into these regions. The first is the available photon energy, Eph, from the active 
nucleus. We use the IRAS 60//m and 100//m flux from Mrk 573 to estimate the nuclear bolometric 
luminosity by assuming that all photons from the central source are absorbed and re-emitted in 
the infrared by surrounding dust, which underestimates the actual luminosity. Each region's Eph 
follows from the fraction of nuclear energy that is intercepted by the arc, the covering fraction c/, 
over its interaction age (iem) assuming that the regions have unity filling factor 



cz 



2 



Eph = 1.5 X 10^^ ^ — j [2.6560 + Sioo] X tem X cf (11) 

For the covering fraction, we assume that the depth of the cloud is approximately the arc thickness, 
9ey The depth could in fact be ~10 times this size, increasing the photon energy by this same 
factor. For our purposes, being conservative with our covering fraction is acceptable, as even with 
this small value, Eph is significantly larger than all other energy sources in the emission regions. 

The second source of energy is relativistic energy, Ey^^i, stored in the radio outflow, which 
provides an estimate of the energy input due to expansion of the radio lobes into the circumnuclear 
gas. -Ej-e/ depends on the minimum magnetic field strength, Bmin, a lower limit found by assuming 
equipartition between the relativistic particles and the magnetic field 

,/ Su a(l + z)^+°^(^)°'-Xo5(a^)\2/7 , ^ 

Brmn « 2.93 X lO-^f— ^^A^_^ii^) (12) 

^"rx"ry J rel^ryy jj^ ) ' 

where X, = 

q — ttr 

Where is the radio spectral index, and frej is the filling factor . The spectral indexes are observed 



to be ar ~ —0.85 for NWl and —0.5 for SEl (jFalcke et al.lll998l ). For the projected size of the radio 



outflow working surface we adopt 6rx and 6ry to be ~0'.'l. We adopt a«2, where a represents the 
contribut ion of relativistic en ergy density by ions in proportion to the contribution by electrons, 
following I Whittle et al.l ([20091) and assume unity filling factor for the relativis tic material (.U pi = 1). 



Following our calculation of B^in, the available relativistic energy, following [Whittle et al.l (|2009l ) 
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is thus ^ 

E,,i ^ 1.6 X 10^6 Q^j2jfj(^^^ Bl,^ (13) 

The estimates for the NWl and SEl regions are hsted in TableSl 

Our calculations of the energies of the different components in the emission line regions and 
the potential sources indicates that Ep}^ dominates, being ~1000 times larger than Ej.(,i and about 
10 times larger than the total energy in the emission regions, -Eem- A complete discussion of the 
different energy components follows in sec. l4.3i 

The energetics of the emission line region and its possible radiative and relativistic sources 
are important for understanding the host-AGN interaction, and how it manifests itself. Similarly, 
examination of the different pressures can reveal whether the relativistic pressure from the radio 
lobe expansion or the radiative pres s ure fr om the central photon source are responsible for the arc 



structures. Following IWhittle et al.l (|2009l ). we calculate the emission-line region thermal pressure 
(Pern)) the relativistic pressure {Prei), and the radiation pressure {Prad), expressed in dynes cm~^: 

Pem ~ 1.4 X 10-9 nem,3TemA (14) 
Prel « OmiBl,J-f (15) 

Prad ~ 2.1 X 10-1^2-6560 + 5ioo]Ae^2 (15) 

where Ag is the distance from the region to the active nucleus. The thermal pressure, Pem, is based 
on the temperatures and densities determined from our spectrophotometric analysis. The relativis- 
tic pressure, Prei, which results from radio lobe expansion is based on the calculated magnetic field 
Bmin- Lastly the radiation pressure from the central photons, Prad, is based on the FIR values, 
similar to our calculation of Ep^. We find Prei-Pem^^5:l, with Prad contributing negligibly. In the 
next section we describe our methods for quantifying the properties of the outflow itself, rather 
than its effect on the emitting material. We combine these two analyses in section lTSl and use the 
results to determine how much influence the nuclear outflow of Mrk573 has on its host. 



4.2. Effect of an Outflow: Velocities and Dynamical Pressures 

Beyond analyzing the energies and pressures of the emission line regions, we can quantify 
parameters of the outflow itself, in particular the velocity and dynamical pressure. This allows us 
to determine the outflow's strength directly, rather than conjecturing it based upon the behavior 
of surrounding material. As the outflow does not appear to be affecting SE2, we only determine 
properties of it with respect to NWl and SEl. 

The strength of the outflow is dependent upon its composition. A relativistic outflow, in the 
form of a jet, with speed Vo approximately equal to the speed of light, is able to both bulldoze 
and shock-ionize material. For an outflow with both thermal and relativistic material, the therma l 



component literally weighs down the outflow, slowing it down to ~ 10'^ km s ^ ( Whittle et al.ll2009l ). 
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Thus, outflow speed reflects the amount of thermal and relativistic material present. [Whittle et al 
(j2009l ) have suggested that the outflow of Mrk78 is a heavy thermal jet, making it too weak to 
have significant influence over its surroundings. Through our quantitative analysis of the outflow 
in Mrk573, we can determine whether or not it is also heavy and weak. 

The approximate speed of the outflow is estimated from the energy and momentum values. 
An outflow with both material components can contribute to the kinetic and relativistic energy of 
the emission line regions; it can also affect the momentum of the emitt ing material, Gpm,. T hus, 



we use these values to determine Vq, the speed of the outflow. Following [Whittle et al.l (|2009l ): 



Vo 



Ekin + E, 



rel 



1 + 



Rkir, 



(17) 



where Rkin is the ratio of kinetic energy to relativistic energy. If we assume the outflow primarily 
consists of heavy thermal material, making Rkin ~ 1 , we calculate velocities of ~ 300 — 2000 km 
s~^(see TableHI, which is greater than the observed emission-line regions' speeds by factors of ~10- 
100. This supports the suggestion that the outflow in Mrk573 is primarily made up of thermal 
material. A primarily relativistic outflow would have velocities vastly more than ~10-100 times 
the emission- line velocities, and we do not calculate speeds of this extent. 

The t hermal component o f the outflow drives dynamical pressure, Po,dym into its surroundings. 
Following [whittle et al.1 



Po,dyn — 



(18) 



where Aq is the area of the outflow, which we measure from the radio maps of lFalcke et al.l (|l998l ) 
to be ~ 0.60 arcsecond^, and we assume that Ho, the force of the outflow, is ~ Hemj the force in the 
emission line regions, calculated using the momentum and crossing times. We can then compare 
the ram pressure of the outflow to the pressures measured for the emission regions and evacuated 
bubbles, finding that it is significantly less than Pem and Prei- Following this, we use the dynamical 
pressure calculation to determine the Mach number of the outflow. 



^Po,dyn 
^Prel 



(19) 



As listed in TablelH the estimated Mach numbers, though uncertain, are very small, on the order 
of 0.01, confirming that t he outflow is slow an d thus, likely dominated by thermal material, similar 
to the outflow in Mrk78 (jwhittle et ahlbood ). 



A final aspect of our outflow analysis is to determine the amount of mass the outflow transports. 
This value quantifies the influence of the outflow on its surroundings. The mass flux is: 



no 

Vo 



(20) 



and is ~ 10 Mq s based on our estimates above. 
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4.3. Discussion: Outflow Strength 

Our spectrophotometric analysis suggested that the outflow in Mrk573 is unable to bulldoze 
material into the emitting arcs. Through analysis of the physical properties of the emission regions 
and the outflow, we can confirm this picture. 

As shown in section HTT] the relative proportions of each of the energy components is 
E^h ■■ E.rr,, : Er,i I Ekir,. : Eth «10'^ : 3000 : 10 : 1 : 0.01 (see Tablei]), similar to what was found in 



Mrk78 by [Whittle et al.l (|2009l ). The photon energy Ep/j clearly dominates. The emission-line 
region energy Eg^ is ~25% of Ep/j, confirming that central photoionization alone is sufficient to be 
the ionizing source of the emission regions. The relativistic and kinetic energies, while both sub- 
stantially smaller than the emission-line region energy, are comparable to each other. This implies 
that the slow expansion of the radio source into the surrounding material is able to drive the small 
amount of kinematic activity in the emission-line regions. Also note that SE2, where there is little 
evidence of interaction between the material and an outflow, has very similar energy magnitudes 
to SEl and NWl. Despite the fact that NWl and SEl are in the direct path of the outflow, the 
regions appear quite similar to non-affected material, implying that an outflow has little influence 
on the emission-line regions' properties. 

We further investigate the kinematics by looking at the different pressures in the regions. The 
radio outflow pressure Prei is slightly larger than the ionized gas pressure, Pem, indicating that 
the small kinematic motions in the emission regions are caused by the gentle expansion due to the 
radio outflows. Our previous analysis of the densities and temperatures of the emission regions and 
the evacuated bubbles shows that the bubbles have comparable or less thermal pressure than the 
arcs over a wide range of temperatures. This supports the conclusion that any expansion of the 
outflows is gentle. The radiative pressure Prad is markedly less than the emission-line region thermal 
pressure, and therefore makes no significant contribution to the emission-line region kinematics we 
observe. We also investigate the dynamic pressure contribution from the thermal component of 



the o utflow, which has been suggested as a potentially significant pressure source (jBicknell et al 



19981 ). However for Mrk 573 the dynamic pressure is less than the relativistic pressure. The dynamic 
pressure would have to be significantly larger than the relativistic pressure or comparable to Pem 
to be able to push around the host ISM significantly. In this case, however, it appears that the 
expansion is caused primarily by the radio lobes, with the mechanical contribution of the outflow 
being supplementary at best, and likely negligible. 

Our calculation of the amount of material transported by the outflow also indicates that it 
has little influence over the host. The outflow can transport 1O~'^M0 yr~^ of material. Over the 
calculated lifetimes of the outflows, they have carried at most a tenth of the emission mass. Strong 



jets have been measured carrying as much as ~O.6M0yr (jBicknell et al.l Il998l ). significantly 
more than we calculate here, once again pointing towards the general impotence of the outflow in 
Mrk 573. 



The energy analysis shows that radio outflow contributes primarily to the kinetic energy of 
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the ionized region but not to its ionization, while the pressure analysis shows that the outflow 
does not strongly expand into the surrounding material. While these estimates have significant 
uncertainties associated with them, the differences between the estimated energies, pressures, etc. 
are large enough that we can clarify the relationship between the outflow and host despite this 
uncertainty. Taken together, the data suggest that the host galaxy of Mrk573 is photoionized 
rather than shock-ionized, and that the radio outflow gently shapes the regions but does not 



otherw ise contribute to its ionization or heating significantly. This is similar to what IWhittle et al 
(|20o3) found for Mrk 78. 



The extensive evidence that the emission regions of Mrk 573 are pre-existing structures pho- 
toionized by a nuclear source also support that the outflow itself is weak and likely thermal. An 
outflow made up of primarily thermal material is slower and therefore less a ble to bulldoze and 
strongly shock surrounding material. Earlier work by iBicknell et al.l (jl998l ) on the outflow in 
NGC 1068 has found similar results, suggesting that the radio jets in Seyferts differ from those in 
radio galaxies and quasars primarily in that they are dominated by thermal gas instead of light 
relativistic material. Our analysis of the outflow speeds confirm that Mrk 573 has a slow outflow, 
orders of magnitude slower than one dominated by relativistic material. This suggests that the 
outflow of Mrk 573 is similar to the weak heav y jet that is incapab le of shocking material and being 
the primary ionization mechanism in Mrk 78 (jWhittle et al.ll2009l ). 



Combining the different facets of this quantitative analysis presents a consistent picture of a 
slow, heavy outflow capable of gently shaping the surrounding host ISM into the arcs we see, but 
not fast enough to substantially heat, displace, or ionize that material. 



5. Summary 

We have examined the influence of the AGN in Mrk 573 on its host galaxy using high angular 
resolution spectrophotometry from HST STIS. These spectra provide both kinematic and spec- 
trophotometric measurements from which we derive the kinematics and physical state of the gas 
in the circumnuclear regions associated with an interaction between a nuclear radio outflow and 
spiral dust arms in the host galaxy. Arguments from emission-line diagnostics and thermodynam- 
ics (energy and pressure in the various components) lead us to conclude that the outflow from the 
active nucleus in Markarian 573 does not strongly influence the surrounding ISM other than gently 
sculpting the material into arcs. The dominant source of the heating and ionization in the gas can 
be entirely explained by photoionization by energetic photons from the active nucleus proper. This 
is similar to what has been seen in other, well-studied nearby AGN outflows with sufficient data 
to make similar energetics arguments, particularly the work of Whittle et al. on Mrk 78. In many 
ways, Mrk 573 and Mrk 78 are very similar energetically and morphologically: they have hot, heavy 
outflows sculpting extended emission-line regions lit-up by ionizing photons from the nucleus. 

All of the outflow and photoionization energetics taken together suggest that while there is 
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some feedback on the AGN's host, it is relatively gentle and insufficient to unbind the host ISM 
or otherwise shutdown star formation by many orders of magnitude. The outflows expand gently 
into the host and have transported at most a tenth of the mass in the emission regions over their 
lifetimes. The star formation regions in the inner few kil oparsecs of Mrk573, seen as wound spiral 
arms outlined by H II regions (jPogge de Robertislfiggi 'l. appear to be undisturbed by the outflow. 
This has important consequences for AGN feedback models in galaxies. Specifically, this shows that 
outflows such as the one in Mrk 573 are not strong enough to remove material from the host galaxy; 
in fac t, analyses of low luminosity AGN indicate tha t AGN feedback has little influence over the 
host (jKrongold et ahlboOTl : Iwhittle et aPbood . bood . bood ). As feedback is frequently invoked to 
explain different aspects of host- AGN evolution, new models must take into account that feedback, 
at least for AGN comparable to local Seyferts, has a significantly smaller effect on the host than 
previously thought. 
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Fig. 1. — Contrast-enhanced structure r nap of the centra. 1 10" o f Mrk573, constructed from an 
archival WFPC2 F606W filter image by iPogge &: Martinil (j2002l ). Emission- Hne regions appear 
bright and dust lanes appear dark. The STIS sht used for the long-sHt spectra shown in Figure [2] 
is superimposed (the actual long slit extends 54", but we sh ow only the inner 1 0"). The main 
emission-line regions are labeled with the names adopted from lFerruit et al.l (jl999l ). The image is 
oriented North up, East to the left, with scale bars indicating the angular and approximate linear 
scales in arcseconds and kiloparsecs, respectively. 



- 24 - 



G430L 

■ ... .i 


— NW1 


. • . . • — SE1 

^ - II III II 1 ■ ' — SE2 
-M- [NeV] [Oil] [Nelll] Hy-i-[0 III] Hell Hp [0 III] 


G750M t t 

i i 4 


t t 


- •» ^ ^ 

1 II • ? ? 

SE [0 1] [0 1] [FeX] Ha + [Nil] 


[Sll] 



Fig. 2.— Long-slit 2D STIS spectra of Mrk573 taken with the G430L (top) and G750M {bottom) 
gratings. The central 10" of the STIS slit is shown, with NW at the top and SE at the bottom 
of each panel. The G430L spectrum extends from A2900A to A5700A, while the G750M spectrum 
extends from A6300A to A6850A. Prominent emission lines are labeled. 

Table 1. Line Flux Measurements 



Feature 


NW 1 


NUC 


SE 1 


SE 2 


Aperture Area (arcsec^) 


0.051 


0.061 


0.041 


0.041 


[Ne V]A3426A 


5.59e- 


17 


3.40e- 


15 


4.04e- 


17 


5.16e-17 


[0 II]A3727A 


S.lle- 


16 


1.51e- 


15 


7.09e- 


16 


1.28e-16 


[Ne III]A3868A 


1.17e- 


16 


1.73e- 


15 


2.42e- 


16 


5.75e-17 


[Ne III]A3968A 


3.76e- 


17 


7.80e- 


16 


1.37e- 


16 


3e-17 


H7 A4340 A 


5.07e- 


17 


8.48e- 


16 


1.03e- 


16 


3.98e-17 


[0 III]A4363A 


l.lOe- 


17 


4.22e- 


16 


3.47e- 


17 


5.90e-18 


He IIA4686A 


2.87e- 


17 


8.96e- 


16 


5.33e- 


17 


1.63e-17 


H/3 A4861A 


1.08e- 


16 


1.85e- 


15 


2.53e- 


16 


7.24e-17 


[0 III]A4959A 


4.70e- 


16 


7.94e- 


15 


9.58e- 


16 


3.10e-16 


[0 III]A5007A 


1.29e- 


15 


2.37e- 


14 


2.80e- 


15 


8.98e-16 


[0 I]A6300A 


5.32e- 


17 


5.37e- 


16 


1.02e- 


16 


1.81e-17 


Ha A6563A 


4.48e- 


16 


7.14e- 


15 


9.76e- 


16 


2.58e-16 


[N H]A6583A 


4.32e- 


16 


3.57e- 


15 


9.62e- 


16 


1.21e-16 


[S H]A6716A 


1.28e- 


16 


6.29e- 


16 


2.69e- 


16 


4.81e-17 


[S H]A673lA 


1.31e- 


16 


8.64e- 


16 


3.24e- 


16 


4.3e-17 



Note. — The line flux measurements for the four emission-line re- 
gions in Mrk573. The fluxes have the units erg cm~^ and uncer- 
tainties of 20% or less. 
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Fig. 3. — The Ha radial velocity profile {top), FWHM profile (middle), and total line intensity 
profile {bottom) along our slit for Mrk573 as a function of distance from the nucleus. The red 
triangles represent rotating disk component, the black squares correspond to the biconical outflow 
region, and the large blue filled circle represents the nucleus. The centroids of the individual 
emission line regions are marked by the dotted vertical lines. 
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Fig. 4. — Same as Figure El for the [O III] A5007A emission line. The FWHM curve is not shown 
as the resolution of the G430L grating is low dispersion, resulting in broader features. 
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Fig. 5. — Schematic of the inferred outflow geometry viewed perpendicular to the line of sight to- 
wards Earth (indicated by the large arrow) . This view is from the south lookin g towards the nucleus. 
The io nization cone is shown with the 45° opening angle {9c) measured by IWilson &: Tsvetanov 
(j 19941 ) and is tilted by angle 18° (0) with respect to the plane of the sky, as estimated from our 
kinematics (see section IXT]) . In this geometry, the NW outflow is emerging out of the plane of the 
sky towards us, while the SE flow is going into the sky. 



Table 2. Density and Temperature Regional Analysis 



Feature NW 1 NUC SE 1 SE 2 



HP 

AA4959+5007 

4363 obs 
AA4959+5007 

4363 corr 
A6716 
A6731 
Uf. (cm^^) 

Te (K) 



0.471 ±0.13 
160.42±54 
160.59±54 

0.979±0.27 
789±78 

10814±1081 



0.458±0.13 
74.96±25 
77.17±26 

0.728±0.20 

2437±243 
14454±1445 



0.407±0.11 
108.36±36 
125.53±42 
0.828±0.23 
1489±148 
12445±1244 



0.549±0.15 
204.90±69 

0.981±0.27 

762±76 
10046±1004 



Note. — The line ratios and derived densities and temperatures for the 
four extraction regions. The Balmer decrement for SE 2 was non-physical for 
the specified temperature, so we assume that the reddening is negligible. The 
uncertainties for the density and temperature are 10%. 
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Fig. 6. — Extracted spectra of emission-line regions in Mrk 573. The left panels contain the G430L 
spectra and are labeled with the associated region from Figure [TJ The right panels contains the 
G750M spectra. Each extraction region was chosen to include all of the emission of each region; the 
angular extent of these regions is listed in Tabled! All of the regions have [O III] AA4363, 4959, and 
5007A and [S II]AA6716 and 6731A emission features, which were used to determine the density 
and temperature. The nucleus includes additional high-ionization features, such as [Fe X] A6375A 
which indicates the presence of close-in material at higher temperatures. 
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Fig. 7. — Optical ionization diagnostic diagrams. These diagnostics are the best at distinguishing 
between shocked and photoionized regions. Inconclusive ionization diagnostics are not included, 
but are available by request. The dashed line is U1.4, photoionization of an optically thick medium 
by power-law emission. The crosses mark log Uia = —1.0, —1.5, —2.0, —2.5 and —3.5. The dot-dash 
line represents \Jdust, power-law photoionized d usty gas. Crosses mark log Udust = —3.0, —2.0, —1.0, 
and 0.0. The solid lines represent shocks from lAllen et al.l (|2008l ). The thin line is a photoionizing 
shock with B/-y/n ~ 0/iGcm~^/^. The thicker solid line connected to the thin is a photoionizing 
shock with B/y^ = 10 fiG cm~^/^. These two shocks are modeled for velocities of 200-1000 km s~^. 
The crosses mark every additional lOOkms"^. These models are shocks and precursors and assume 
a pre-shock density of 1 cm~'^ and solar abundance. The points represent each of the four studied 
regions of Markarian 573. Blue represents the nucleus, green is the NWl region, and red and black 
are SEl and SE2, respectively. 
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Table 3. Emission Measure Analysis 





Temperature (K) 






I(H/3) (erg s^^ cm^^ arcsec^^) 




4.1e-16 


2.4e-16 


Depth (arcsec) 




1.6 


1.6 


Depth (parsec) 




576 


576 


Emission Measure (cm^^ pc) 


15000 


826 


483 


(rie) (cm"3) 


15000 


1.2 


0.9 


Emission Measure (cm~^ pc) 


10^ 


36210 


21197 


(ne) (cm-3) 


10^ 


7.9 


6.1 



Note. — Temperature and density calculations from the emission mea- 
sure for the bubble regions as described in section 13.2.21 B nw refers to 
the bubble between NWl and the nucleus and Bse between the nucleus 
and SE 1. 
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Table 4. Quantitative Analysis of the Outflow Arcs 



Quantity 


Symbol (umtsj 


i\ VV 1 


In U (.^ 


ash i 




Region Angular Size (x) 


9ex (arcsec) 


2.14±0.1 


0.67±0.1 


2.55±0.1 


3.72±0.1 


Region Angular Size (y) 


O^y (arcsec) 


0.31±0.1 


0.31±0.1 


0.26±0.1 


0.31±0.1 


Distance from Nucleus 


Ae (arcsec) 


1.81±0.1 




1.76±0.1 


2.63±0.1 


Electron Temperature 


Tp (K) 


10814±1081 


14454±1445 


12445±1244 


10046±1004 


Electron Density 


rig (cm~^) 


789±78 


2437±243 


1489±148 


762±76 


H/3 Emission Flux 


loa( Fmr) fere; s^^ cm~^) 


-14.8±0.1 


-14.0±0.1 


-14.4±0.1 


-14.9±0.1 


[0 III]A5007 Emission Flux 


Ioq( F^nm) fere: cm~^) 


-13.7±0.1 


-12.9±0.1 


-13.4±0.1 


-13.9±0.1 


Emission-line Width 


W^m Hrv fkm/s) 


159±20 


217±20 


214±20 


108±20 


Emission-line Velocity 


Vim Hrv fkm/s) 

' GTfl jii CK V / / 


110±20 




-83±20 


-62±20 


Emission Mass 


Mp„ (10^Mn)±22% 

J. r J. fJ_jyi y-M-\j J. r J. y _i_ i_i i_i / vj 


0.009 


0.018 


0.012 


0.007 


Emission Region Age 


tem (Myr) 


2.24±0.65 




3.53±1.31 


6.90±2.50 


Emission Luminosity 


looiLprn) fereVs) 


41.1±0.1 


41.9±0.1 


41.4±0.1 


40.9±0.1 


Emission Region Energy 


Ioq( E^m) ferg) 


54.97±0T5 




55.46±0.18 


55.26±0.18 


Emission Region Momentum 


looiG^rr,^ fgm cm s~^) 


44.81±0.15 




44.81±0.18 


44.45±0.18 


Emission Region Force 


/oc/fn^m) fdvne) 

J \ Gill } V -/ / 


30.96±0.22 




30.76±0.26 


30.11±0.26 


Thermal Energy 


loq(Eth) (erg) 


49.60±0.10 


50.03±0.10 


49.79±0.10 


49.46±0.10 


Transverse Kinetic Energy 


loqiEkjn f) fere;) 


52.06±0.15 




51.94±0.18 


51.45±0.18 


Internal Kinetic Energy 


log{Ekin,i) (erg) 


50.61±0.10 




51.00±0.10 


50.17±0.10 


Total Mechanical Energy 


log{Emec) (erg) 


52.07±0.18 




51.99±0.20 


51.47±0.20 


Photon Energy 


log{Eph) (erg) 


55.6±0.1 




55.9±0.2 


56.1±0.2 


Relativistic Energy 


log{Erei) (erg) 


«52±0.17 




«53±0.17 




Emission Region Pressure 


log{Pe,n) (dyne cm^^-j 


-8.92±0.04 


-8.31±0.04 


-8.59±0.04 


-8.97±0.04 


Radiation Pressure 


logiPrad) (dyne cm^^-j 


-10.5±0.06 




-10.5±0.06 


-10.9±0.04 


Relativistic Pressure 


log{Prei) (dyne cm"^) 


-8.0±0.14 




-7.6±0.14 




Outflow Ram Pressure 


log{Po.dyn) (dyne cm^^) 


-10.9±.22 




-11.1±.26 




Outflow Speed 


K (km/s) 


340±228 




1700±1207 




Outflow Mach Number 




0.03±0.5 




0.01±0.5 




Magnetic Field Strength 


log{Bjnin) (Gauss) 


-3.25±0.07 




-3.05±0.07 




IRAS eO^fm Flux 


•S'eo (Jy) 




1.27±0.15 






IRAS lOO^im Flux 


S'loo (Jy) 




1.26±0.16 







Note. — A quantitative analysis of Mrk573 calculated based upon the formulations of I Whittle et al. 
quantities are derived in section [4l Uncertainties are calculated using standard error propagation. 
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